SELF-CONCORDANCE IS NP-HARD 



LEK-HENG LIM 



Abstract. We give an elementary proof of a somewhat curious result, namely, that deciding 
whether a convex function is self-concordant is in general an intractable problem. 



1. Introduction 

Nesterov and Nemirovskii [20] famously showed that the optimal solution of a conic programming 
problem can be computed to e-accuracy in polynomial time if the cone has a self-concordant 
barrier function whose gradient and Hessian are both computable in polynomial time. Their work 
established self-concordance as a singularly important notion in modern optimization theory. 

We show in this article that deciding whether a convex function is self-concordant at a point 
is nonetheless an NP-hard problem. In fact we will prove that deciding the self-concordance of 
a convex function defined locally by a cubic polynomial (which cannot be convex on all of R"), 
arguably the simplest non-trivial instance, is already an NP-hard problem. This is not so surprising 
given that deciding the convexity of a quartic polynomial, arguably the simplest non-trivial instance 
of deciding convexity, is also NP-hard — a much harder recent result [T]. In addition to the NP- 
hardness of self-concordance, we will see that there is no fully polynomial time approximation 
scheme for the optimal self-concordant parameter and that deciding second-order self-concordance 
[14) of a quartic polynomial is also an NP-hard problem. 

These hardness results are intended only to add to our understanding of self-concordance. They 
do not in anyway detract from the usefulness of the notion since in practice self-concordant barriers 
are constructed at the outset to have the requisite property [20l Chapter 5]. It is unlikely that one 
would ever need to generate random functions and then test them for self-concordance. 

We deduce the NP-hardness of self-concordance using a result of Nesterov himself, namely, 
minimizing a cubic form over a sphere is in general an NP-hard problem. Nesterov's result, which 
appears in an unpublished manuscript [19], contains some minor errors that have unfortunately 
been widely reproduced. We take the opportunity here to correct them and also to prove a variant 
of Nesterov's result directly from Motzkin-Straus Theorem [18| . 



2. Self-concordance in terms of tensors 

Let Q C M" be open and / : 17 — ?• M be in C"^(0), i.e., has continuous partials up to at least 
order d. Recall that the dth. order derivative of f ai x £ Q,, denoted V°'/(x), is a tensor of order d 
[15) . To be more precise, this simply means that V^fix) is a multilinear functional on Tx{^), the 
tangent space of at x, that is, 

V'^fix):T,{n)x---xT,{n)^R, 

V 

d copies 
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where 



V'^/(x)(/ii, . ..,ahi + ph'i, ...,hn) = aV^fix^hi, . . . ,hi, . . . ,K) 



+ /3V'^/(x)(/ii,...,/i^,...,/i„) fori = l, 



, n. 



With respect to the standard basis . . . , of Tj(ri), we may identify T^{n) ^ M" and V^fix) 
may be regarded as a 'd-dimensional matrix' (d-hypermatrix for short), 

V'fix) = [a.,...,]r„....,=i e M-x-x". 

Indeed, we must have 

^ ^'/(^) 
dxi-^ ' ' ' dxi^ 

and since / G C"^(r2), we get that ai^...i^ = «j^(i)---V(d) ) V'^/(a;) is a symmetric d-hypermatrix. 
Every symmetric d-hypermatrix A = [aii-.-i^]"^ G i^"-x-"X" defines a homogeneous polynomial 

of degree d, denoted 



A{h,...,h) := ^ ai^...i^hi^---hi^£M\hi,... ,hn_ 



d- 



=1 



Readers should be able to infer from the above discussion that d-hypermatrices are coordinate 
representations of d-tensors, just as matrices are coordinate representations of linear operators and 
bilinear forms (both are 2-tensors). We refer the reader to pTJ for more information. 

The usual definition of self-concordance requires that / S C'^{Q) and in which case it is given by 
a condition involving the matrix V^/(x) G M"^" and the 3-hypermatrix V^/(x) S ^nxnxn^ 

Definition 2.1 (Nesterov-Nemirovskii) . Let 0, C R" is a convex open set. Then / : — )• M is said 
to be self- concordant with parameter a > at x £ if 

(2.1) V^fix)ih,h)>0 
and 

(2.2) [V'fix)ih, h, /i)]2 < 4a[VV(x)(/i, h)f 

for all h G M"; / is self- concordant on 0, if (j2.ip and (j2.2p hold for all x G 0. The set of 
self- concordant functions on 17 with parameter a is denoted by Sa{^)- 



By (12. a function self-concordant on $7 is necessarily convex on Q. A minor deviation from 
|20j is that cr above is really the reciprocal of the self-concordance parameter as defined in |20^ 
Definition 2.1.1]. Our hardness results would be independent of the choice of a. Note that 

i,j=l ■' i,j,k=l ■' 

So for a fixed x G fi, V'^f{x){h, h) is a quadratic form in h and V^f{x){h, h, h) is a cubic form in 
h. It is well-known that deciding (j2.ip at any fixed x is a polynomial-time problem (but not so for 
deciding it over all x £ see [Ij). Hence given a cr > 0, deciding self-concordance at x essentially 
boils down to (j2.2p : Is the square of a given cubic form globally bounded above by the cube of a 
given quadratic form? We shall see in the next sections that this decision problem is NP-hard. 

While we shall think of V^/(x) as a matrix and V'^f{x) as a hypermatrix thoughout this article, 
we nonetheless wish to highlight that condition ()2.2p is really a condition on V^/(a;) regarded as a 
2-tensor and V'^/(x) regarded as a 3-tensor; that is, (12. 2p is independent of the choice of coordinates, 
a property that follows from the affine invariance of self-concordance [20', Proposition 2.1.1]. Self- 
concordance on Q is then a global condition about the tensor fields V^/ and V^/. 
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3. Maximizing a cubic form over a sphere 



We will include a proof that the clique and stability numbers of a graph with n vertices and m 
edges may be expressed as the maximal values of cubic forms (in n+m variables) over the unit sphere 
gn+m-i_ "Yhis, or at least the stability number version, is known but the reference |19l Theorem 4] 
usually cited contains some slight errors that have been reproduced several times elsewher^ll- We 
take the opportunity to provide a corrected version below. Our proof follows Motzkin-Straus 
Theorem [18j and the similar result of Nesterov [191 Theorem 4] for stability number. 

Let G = {V, E) be an undirected graph with n vertices and m edges. We shall require that 
E ^ throughout, so n > 2 and m>\. Recall that S" C 1/ is a clique in G if € E for all 

i,j G S and S C 1/ is stable in G if ^ E for all i,j € S. The clique number and stability 

number of G are respectively: 

oj{G) = maxjl^l : 5 C y is clique}, a{G) = max{|S'| : 5 C y is stable}. 

Motzkin and Straus [18] showed that uj{G) may be expressed in terms of the maximal value of a 
simple quadratic polynomial over the unit simplex. Although not in [18j, it is straightforward to 
see that essentially the same proof also yields a similar expression for a{G). 

Theorem 3.1 (Motzkin-Straus). Let A" = {x € M" : xi + • • • + a;„ = 1, > 0} denote the unit 
simplex in M". Then the clique number oo{G) and stability number a{G) may be determined via 
quadratic optimization over simplices: 

^ Y ^ 

(3.1) 1 7--r = 2 max > XjX,-, 1 -— - = 2 max > 



3- 



u{G) xeA"^{i,j}eE J' a{G) 

Since deciding if a clique of a given size exists is an NP-complete problem [11], an immediate 
consequence is that clique number is NP-hard and by Motzkin-Straus Theorem, so is quadratic 
maximization over a simplex. 

In an unpublished manuscript [19^ Theorem 4], Nesterov used Motzkin-Straus Theorem to obtain 
an alternate expression ()3.3I) for stability number involving the maximal value of a cubic form 
over a sphere. In the following we derive a similar expression (13. 2p for the clique number, which 
yields slightly simpler expressions for our discussions in Sections [5] and [Tj and may perhaps be of 
independent interest. 

Theorem 3.2 (Nesterov). Let G = {y,E) be an undirected graph with n vertices and m edges. Let 
Ei'^^^ = {x G M'^ : ||x|| = 1} denote the unit I'^-sphere in Mf^. The clique number uj{G) and stability 
number a{G) may be determined via cubic optimization over spheres: 



(3.2) 



(3.3) 



1 



a;(G) 
1 

a{G) 



27 

— max 

2 (u,«)eS"+'"-i 

27 

— max 

2 {u,v)&"+"^-'^ 



'{i,j}eE 



UiUjWij 



UiUjWij 



Proof. This follows from Motzkin-Straus Theorem and the equalities 
(3.4) max > x,Xi = max > 



^{i,j}eE 



2 2 

'{i,j}eE * J 



(3.5) 
(3.6) 



max 

27 

— max 

4 (n,i;)e§"+'"-i 



E 



UiUjWij 



^{i,j}(^E 



UiUjWij 



^Foi example, [51 Theorem 3.4]. To see that the expression is incorrect, take a graph with three vertices and one 
edge, the left-hand side gives and the right-hand side gives 1. 
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3.4p comes from substituting Xi = uj, i = 1, . . . ,n. Cauchy-Schwartz yields 

/ N 1/2 / N 1/2 

{i,j}&E * J *J - yZ^{ij}eE ' 3 J \Z^{i,j}(^E J 



and so 



max 



||u|| = ||ty||=l^{ij}e-B 



UiUjWij < max 



(3.7) 



max ( >^ ufu^i I max ( >^ wf^ ] 

u\\=i \^{iJ}eE ' 3 J 11^11=1 \Z^{i,j}eE J 

/ \ 1/2 

max ( >^ ufui I =: a. 

u\\=i \^{i,j}eE ' J J 



Let the maximal value a be attained at u G S" We set Wij = uiUj/a for all € E (note 

that a > if / 0). Observe that 



EUiUjWii = — ufu'^ = a, 

{i,j}eE ^ a'^iiJjeE ' ^ 



and w G S™" ^ since 



2 1 \ ^ 2 2 

Hence equality is attained in p.7p and we have (|3.5p . We deduce (|3.6p from 



^ _IU^, = — > , , = 1. 



max 



: sup 

/3e(o,i) 

: sup 

/3e(o,i) 



{i,i}ei? 



max 



. UiUjWij 



max 



. UiUjWij 



1-/3 ^{ 



||u||2=^^ ll^l|2=i_^^{i,j}e£; 



UiUiW. 



max 



Ui Uj Wij 



max > UiUiWii 
^,||2=i, ||^||2=iZ^{i,i}e£; ^ 

2 

— ^ max > UiUnWii 

3^/3\\u\\ = \\u,\\=l^{i,j}&E ^ 



X sup — (3 

/3e(o,i) 



Note that the maximal value of X^jj WjUjii^jj, whether over x S™ or over §"■+'" can 
always be attained with n > and w > 0, thereby allowing one to take square in (j3.5p and 
p.6p . The same proof works word-for-word for stability number with the replacement of index of 
summation G E' by '{«, j} ^ -E'. □ 



4. Complexity theory for casual users 

In this article, we use complexity classes defined in the standard classical framework: Time 
complexity measured in bits with the Cook-Karp-Levin notions of reducibility [H \TT\ [16] on a 
Turing machine |24j . This is also the most common framework for discussing complexity issues in 
optimization [25]. We briefly recall the intuitive ideas behind the complexity classes used in this 
article for our readers. This is of course not meant to be anywhere near a rigorous treatment; for 
that, see [m [23]. 

A decision problem, i.e., answer is yes or NO, is in NP if one can verify an yes answer in 
polynomial time; it is said to be in co-NP if one can verify a NO answer in polynomial time. A 
decision problem in NP is said to be NP-complete if one can reduce any other problem in NP to 
it. Likewise, a decision problem in co-NP is said to be co-NP-complete if one can reduce any other 
problem in co-NP to it. NP-complete and co-NP-complete problems are believed to be intractable. 
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Example 4.1 (Subset sum problem [^). Let S be a finite set of integers. The problem "Does S 
have a non-empty subset with a zero sum?" is NP-complete — one can check a purported yes 
answer, i.e., a nonempty subset of integers, is indeed a yes answer, i.e., has a zero sum. The 
problem "Does every non-empty subset of S have a nonzero sum?" is co-NP- complete — one can 
check a purported NO answer, i.e., a nonempty subset of integers, is indeed a NO answer, i.e., has 
a nonzero sum. Note that the two problems are logical complements of each other. This is in fact 
another way to define the co-NP class, namely, it comprises problems that are logical complements 
of NP problems. 

A problem is said to be NP-hard 112\ [13] if one can reduce any NP-complete decision problem 
to it. A problem is said to be co-NP-hard if one can reduce any co-NP-complete decision problem 
to it. In other words, if one can solve an NP-hard problem, then one can solve any NP (including 
NP-complete) problems; if one can solve a co-NP-hard problem, then one can solve any co-NP 
(including co-NP-complete) problems. So NP-hard problems are believe to be even harder than 
NP-complete problem and co-NP-hard problems are believed to be even more intractable than 
co-NP-complete problems. An NP-hard problem need not be in NP, i.e., one need not be able 
to check an yes answer in polynomial time; an NP-hard problem that is in NP is by definition 
NP-complete. Likewise a co-NP-hard problem need not be in co-NP, i.e., one need not be able 
to check a NO answer in polynomial time; a co-NP-hard problem that is in co-NP is by definition 
co-NP-complete . 

We have used the term 'reduce' without stating what it meant. There are in fact two different 
notions: The Cook reduction (also known as polynomial-time Turing reduction) and the Karp 
reduction (also known as polynomial-time many-one reduction). For our purpose, all we need to 
know is that Cook reduction is believed to be a stronger notion of reducibility than Karp reduction: 
Under Cook reduction, every decision problem can be reduced to its complement. A consequence is 
that with respect to Cook reducibility, there is no distinction between NP and co-NP, NP-complete 
and co-NP-complete, or NP-hard and co-NP-hard. Under Karp reduction, it is believed (although 
unknown) that these classes are all distinct. 

Given an NP-hard maximization problem, e.g., the ones in Theorems 13.11 and 13.21 a polynomial 
time approximation scheme (ptas) is an algorithm that, for any fixed e > 0, would produce an 
approximate solution within a factor of 1 — e of the maximal value, and has running time polynomial 
in input size and e. 

5. Complexity of deciding self-concordance 

The recent resolution of Shor's conjecture by Ahmadi, Olshevsky, Parrilo, and Tsitsiklis [Ij 
shows that deciding the convexity of a quartic polynomial globally over is NP-hard. So the 
self-concordance of a function that is not a priori known to be convex is NP-hard in a trivial way 
since deciding whether (12. ip holds for all x £ 0, in Definition 12.11 is already an NP-hard problem. 
Our complexity result requires more stringent conditions: (i) Our functions are assumed to be 
convex in and so (12. ip is always satisfied and self-concordance reduces to checking (12. 2p : (ii) We 
show that (j2.2p is already NP-hard to check at a single point x E fi. 

Throughout the following we shall require the inputs to our problems to take rational or finite- 
extension^ of rational values, e.g., A £ Q"^"^",q E Q, to ensure a finite bit-length input. Note 
however that an NP-hard problem may not be in the class NP; so an NP-hard decision problem can 
be posed over the reals, e.g., 'Is there an /i E M" such that [A{h, h, h)]"^ < q[h~^ h]^ holds?' without 
causing any issue since it is not required to have a polynomial-time checkable certificate. 



In this article we only encounter simple quadratic and cubic extensions, cf. (|7.4p and (|5.3|l . Note that for positive 
g G Q, elements of Q(y^) and Q(^/g) may be written as a + b^ and a + b^/g + c(^/g)^ respectively. Therefore they 
may be represented by pairs and triples of rational numbers. 
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We will now formulate a decision problem that will lead us to the requisite NP-hardness of self- 
concordance. Let G = (y, E) be an undirected graph with n vertices and m edges where n > 2 and 
m > 1. Let ^ = {{ik,jk} ■.k = l,...,m]. Define Aq = KkYl^^=i G Q(n+m)x{„+m)x(n+m) 

otherwise. 



^ijk 



Note that Ac is a symmetric hypermatrix, i.e., 

Oiijk — Oiikj — Oijik — (^jki — (^kij — (^kji 

for all i, j,k £ {1, . . . ,n + m}. Let us denote the coordinates of /i G 1 

h= {ui,...,Un,Wi,j,,...,Wi^j^). 

In which case, 



as 



AG{h,h,h) = Y,'^^^ 



y 



Hi Uj Wij , 



and so by Theorem 
(5.1) 



max 



AG{h,h,h) 



max[^c(/.,/.,/.)]^ = -^l--^ 



The CLIQUE problem asks if for a given graph G and a given fc G N, whether G has a clique of 
size k? CLIQUE is well-known to be NP-complete [11]. In other words, deciding if uj{G) > k, or 
equivalently, uj{G) > /c — 1, is an NP-hard problem; and by (|5.ip . so is deciding if 

2 



(5.2) 



max 



AG{h,h,h) 
\\h¥ 



1 



k-1 



Strictly speaking, we have restricted ourselves to the subclass of undirected graphs with at least 
two vertices and one edge — clearly clique is still NP-complete for this slightly smaller class. But 
with this restriction, we may assume that k > 2 and therefore the right-hand side of (|5.2p is always 
defined. 

Hence deciding if there exists an h £ M""*""^ for which 



[AGih,h,h)f > 



27 



1 



1 



k-1 



is an NP-hard problem. Note that q = ^[1 — {k — 1) ^] G Q and so this problem is of the form: 

Problem 5.1. Given a symmetric A G q{"+''")x("'+™)x{"+™) and a positive q £ Q, is it true that 
there exists h G for which [A{h,h,h)]^ > q[h'^hf? 



Let cr G Q, cr > 0, be a self-concordance parameter and let 



1 



1/3 



We follow the notation in Section[2j Let Q be the e-ball i?e(0) = {x G M" : ||x|| < s} where e > is 
to be chosen later. We are interested in deciding self-concordance at a; = of the cubic polynomial 
/ : O ^ M defined by 

f[x) = -x X + Ag{x,x,x) = - y + aijkXiXjXk- 

We have V^/(0) = 7/ where / is the (n -|- m) x (n -|- m) identity matrix. Since 7 > 0, i.e., V^/(x) 
is strictly positive definite in a neighborhood of x = and so there exists some -Be(O) on which / 
is convex — this gives us our choice of e. Also, V'^/(0) = Ag- 
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Hence V'^f{0){h,h) = -^h^h = 7||/i|||, V^f{0){h,h,h) = AG{h,h,h), and / is self-concordant at 
the origin with parameter cr G Q if and only if 

[Acih, h, h)]' < Aa^'ih-^h]' = ^ (l - ^) t^^^]' 
for all h G Note that this problem is of the form: 

Problem 5.2. Given a symmetric A G Q{n+in)x{n+m)x(n+m) ^ positive q £ Q, is it true that 
for every h € R"+"', we have [A{h,h,h)]'^ < qlh'^h]^? 

Mathematically, Problems 15.11 and 15.21 are of course equivalent, being logical complements of 
each other. However they may or may not have the same computational complexity. By our 
discussion in Section [H if our notion of reduction is the Cook reduction, then we may indeed 
conclude that Problems 15.11 and 15.21 are equivalent in terms of computational complexity, i.e., 
deciding self-concordance is NP-hard. However, if our notion of reduction is the Karp reduction, 
then what we may deduce from the NP-hardness of Problem 15. II is that Problem 15.21 is co-NP-hard. 
In either case, our conclusion is that self-concordance is intractable. 

Theorem 5.3. Deciding whether a cubic polynomial is self- concordant at the origin is NP-hard 
under Cook reduction and co-NP-hard under Karp reduction. 

The argument in this section clearly works not just for cubic polynomials but for any / G C^{Q,) 
as long as G r?, V^/(0) = 7/, and V^/(0) = Aq — other derivatives and the remainder term in 
the Taylor expansion of / at x = may be chosen arbitrarily as long as / stays convex in 0. This 
liberty allows one to extend the construction above to functions with other desired properties. For 
instance, we may want an example where = M" and since cubic polynomials cannot be convex 
on the whole of M", we will need a quartic / and therefore need to choose V^/(0) accordingly; or 
we may want an example where / is a barrier function, which is equivalent to / having an epigraph 
{{x,t) G M"+^ : X en, fix) < t} that is closed. One may trivially replace by any point a G M" 
by considering the function fa{x) = f{x — a) on Q = B^^a). 

While we have proved our hardness result for functions on C R", it is easy to extend this to 
any M- vector space, for example, symmetric matrices g"^"- or polynomials M[xi, . . . ,Xn], or even 
Riemannian manifolds with a non-trivial class of geodesically convex functions (i.e., not just the 
constant functions). Since self-concordance at a point is a local property, a choice of coordinate 
patch would transform the problem to one over R"; and by our remark at the end of Section [21 it 
will in fact be independent of our choice of coordinates. 

Deciding self-concordance on the whole of ft is of course at least as hard as deciding self- 
concordance at a point in Q and hence we have the following. 

Corollary 5.4. For any $7 and any a > 0, deciding membership in is NP-hard. 

One may wonder why our conclusion in Theorem 15.31 is stated as NP- and co-NP-hardness as 
opposed to NP- and co-NP-completeness. It may appear that given a NO certificate h G M", it would 
be easy (i.e., requires polynomial time) to decide whether [A{h, h, h)]"^ < q[h~^ h]^ is indeed violated. 
But observe that it is only easy to compute the quantities [A{h,h,h)]'^ , q[h~^h]^, and compare 
their magnitudes when we measure time complexity in units of real operations (i.e., arithmetic and 
ordering in M). Since we measure time complexity in units of bit operations, even if the certificate 
h is in Q", it could well have an exponential number of bits and thus it is not at all clear that we 
may check [A{h, h, h)]'^ < q[h~^ h]^ easily. 

6. Inapproximability of optimal self-concordance parameter 

Let A G Q"^"^" be symmetric and / : — )• M be defined by the cubic polynomial f{x) = 
^x~^x + A{x,x,x). As in Section [SJ O is chosen to be a neighborhood of the origin so that / is 
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convex on Jl. The condition (j2.2p for self-concordance of / at j; = with parameter a > may be 
written as 

(6.1) \Aih,h,h)\<2V^\\h\\l 

for all h G M". This is equivalent to requiring 

, X A(h,h,h) r- 

6.2 max \ ' ' ' < 2^, 

h^o \\h\\l 



as A{—h, —h, —h) = —A{h, h, h) and we may drop the absolute value in ()6.ip . 
Since A £ i^"-x'^x" ig a symmetric 3-hypermatrix, the spectral norm [71[l7] of ^, 

^(/ii,/i2,/i3) A{h,h,h) 
^ 2,2,2:= max II ■ II II II = max — -—p, — . 

hiMM^^ ||'il||2||'i2||2||'i3||2 h^'^ \\h\\2 

For the interested reader, the second equality above follows from Banach's result on the polarization 
constant of Hilbert spaces [21 [22]. Hence the optimal self-concordance parameter of / at a; = 0, i.e., 
the smallest value of a so that ()6.2p holds, is given by 

(6-3) O-opt = ^P||i,2,2- 

The spectral norm of a 3-hypermatrix is NP-hard to approximate to within a certain constant 
factor by [71 Theorem 1.11], which we state here for easy reference. 

Theorem 6.1 (Hillar-Lim) . Let A e Q"xnxn he the input size of A in hits. Then it is 

NP-hard to approximate ||^||2,2,2 to within a factor of 1 — e where 

^ = ^~ i^^ N{N -I)) " 2N[N -I) ^'^{m 

By (j6.3p and Theorem 16.11 (Topt is NP-hard to approximate to within a factor of |(1 — e)^ and 
consequently we have the following inapproximability result. 

Corollary 6.2. There is no polynomial time approximation scheme for determining the optimal 
self- concordance parameter Uopt unless P = NP . 

We refer the reader to [U [9] for more extensive approximability results and approximation algo- 
rithms (that are not ptas). In particular, the results in [9] for quartic polynomials would apply to 
the optimal second-order self-concordance parameter (see the next section). 



7. Complexity of deciding second-order self-concordance 

There is also an interesting notion of second-order self-concordance due to Jarre [14j . This 
requires that / G C"^(il) and is given by a condition involving the matrix V^/(a;) £ R"^" and the 
4-hypermatrix V^fix) £ M"X"X"X". 

Definition 7.1 (Jarre). // C R" is a convex open set, then f : 0, ^ M is said to be self- 
concordant of order two with parameter t > at x £ if 

(7.1) V''fix)ih,h)>0 
and 

(7.2) V^f{x){h, h, h, h) < 6r [V'^f{x){h, h)] ^ 

for all h G R"; / is self- concordant of order two on O if (I7.ip and (17. 2p hold for all x £VL. 
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Note that 

V^f(x)(h, h, h, h) = ^ ^y^^\ hihihkhi, 

is a quartic polynomial in h for any fixed x G Q. 

We follow the same argument in Section [5] to show that deciding (j7.2p is NP-hard. This time 
the result would be deduced from Motzkin-Strass Theorem except that for better parallelism with 
Section [5l we will use the quartic-maximization-over-sphere form (j3.4p instead of the quadratic- 
maximization-over-simplex form (j3.ip . 

Given a graph G = {V,E) with n vertices and m edges where n > 2 and m > 1, we define 
Ag G M"X"X"X" by 



ijkl 



1 i = k, j = I, and {i,j} E -E, 
otherwise. 



So yl = [aijki]'ij fc G Q"X"X"'X»^ ig a symmetric 4-hypermatrix. Now observe that, as in ([3 



max jdcf/i, /i, /i, /i) = max > 



= ^ ( 1 



1 / 1 



by Motzkin-Strass Theorem. As in Section [5l given an integer k >2, deciding if a fc-clique exists 
in G is equivalent to deciding if uj{G) > k — 1. Hence deciding if there exists h G with 

(7.3) Acih, h, h,h)>^(^l- [h'hf 

is NP-hard. 

Given a second-order self-concordance parameter r € Q, r > 0, let 



(7.4) 7- 



1/2 



12r V A; - 1, 

We may now define / : $7 — )• M accordingly as the quartic polynomial 

/(x) = -X X + Ag{x, X, X, x) = - 2^ + 2^ . ^. ^ aijkiXiXjXkXi. 

Hence V^f{0){h,h) = jh^ h = 7||/i||| and V^f{0){h,h,h,h) = AG{h,h,h,h). Again we choose 
to be a neighborhood of the origin so that / is convex on 0, as we did in Section [5j So the function 
/ is second-order self-concordant at x = with parameter r if and only if 

(7.5) [AGih,h,h,h) < GT^^ih-^hf = ^ (l - [h'^hf 

is satisfied for all h E M". As in Section [H we observe that the problem of deciding if there exists 
an h £ M" satisfying (j7.3p and the problem of deciding if (j7.5p is satisfied for all h G M" are logical 
complements. Since the former is NP-hard, we arrive at the following conclusion: 

Theorem 7.2. Deciding if a quartic polynomial is second-order self- concordant at the origin is 
NP-hard under Cook reduction and co-NP-hard under Karp reduction. 



It has recently been shown that deciding various seemingly innocuous properties of quartic poly- 
nomials m [To] all fall into the NP-hard category. Theorem 17.21 provides yet another such example. 
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8. Conclusion 



As we have mentioned in the introduction, the hardness results here are intended to shed Ught on 
the properties of self-concordance. They do not in anyway invalidate the usefulness of the notion in 
practice since there are basic principles that one may use to construct self-concordant (and second- 
order self-concordant) functions for use as barriers in cone programming — see "How to construct 
self-concordant barriers" in \20\ Chapter 5] for an extensive discussion or "Self-concordant calculus" 
in [31 Section 9.6] for a summary. 

Self-concordance and second-order self-concordance are conditions involving high-order tensors 
(orders 3 and 4 respectively), which is a topic of great interest to the author. In particular, their 
NP-hardness serves as yet reminder of the complexity of tensor problems [7]. 

9. Acknowledgement 
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